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We introduce a mesoscale technique for simulating the structure and rheology of block copolymer 
melts and blends in hydrodynamic flows. The technique couples dynamic self consistent field theory 
(DSCFT) with continuum hydrodynamics and flow penalization to simulate polymeric fluid flows 
in channels of arbitrary geometry. We demonstrate the method by studying phase separation of an 
ABC triblock copolymer melt in a sub-micron channel with neutral wall wetting conditions. We 
find that surface wetting effects and shear effects compete, producing wall-perpendicular lamellae 
in the absence of flow, and wall-parallel lamellae in cases where the shear rate exceeds some critical 
Weissenberg number. 



Predicting the non-equilibrium dynamics of multi- 
component inhomogeneous polymer formulations is im- 
portant for development and improvement of paints, ad- 
hesives, cosmetics, and processed foods. Many such sys- 
tems have fluid structures with length scales spanning 
nanometers to millimeters, and stabilization of the finest 
scale structures is often assisted by the presence of block 
or graft copolymers 1]. Computational tools which pre- 
dict equilibrium properties of polymeric materials [3, |^ 
have been developed as a cost effective alternative to 
trial and error experimentation. However, bulk proper- 
ties such as strength, optical clarity, and elasticity often 
depend in detail on the size and distribution of inhomo- 
geneities in the system, which are determined both by 
the composition of the fluid and the manner in which it 
was processed. 

We have developed a mesoscale technique aimed at ex- 
amining the effects of hydrodynamic transport on poly- 
mer self assembly and viscoelasticity in industrial pro- 
cessing flows. The scheme couples dynamic self consis- 
tent field theory (DSCFT) 0, H H with a multiple- 
fluid Navier- Stokes system and a set of constitutive equa- 
tions which describe viscoelastic stresses in a polymeric 
fluid. Both the thermodynamic and hydrodynamic por- 
tions of the model employ rigid wall fields which represent 
a slightly porous material that may be shaped to form 
channels with arbitrarily complex geometries. Each rigid 
object in the system may be fixed or moving, simulating 
the interaction of machine elements with polymers in the 
melt. We refer to the scheme as hydrodynamic self con- 
sistent field theory (HSCFT) to underscore its emphasis 
on hydrodynamic transport. 

In contrast to "phase field" techniques ^ 0, , this 
method simulates the assembly of multiblock copolymer 
mesophases where the polymeric nature of the chains is 
explicitly taken into account. It also differs from SCFT 
and DSCFT techniques in its ability to model hydrody- 
namic transport in complex channels with moving ma- 



chinery. Generally speaking, SCFT is capable of describ- 
ing equilibrium morphologies and micro-phases bound- 
aries. DSCFT is capable of describing non-equilibrium 
systems in which hydrodynamic transport may be ne- 
glected, including phase separating melts and systems 
subjected to simple shear fields. HSCFT, in contrast, is 
appropriate for the description of non-equilibrium sys- 
tems in which hydrodynamic effects play an important 
role including dynamic polymer nanocomposites, mi- 
crofluidic channel flows, and micro-scale rheometry of 
material properties. 

In this letter, we introduce a multi- fluid generalization 
of the "two-fluid" model of Doi and Onuki followed 
by a summary of the thermodynamic equations used to 
model block-copolymer self assembly in the presence of 
rigid boundary surfaces. We illustrate the method by 
examining the effects of hydrodynamic transport on a 
phase separating triblock copolymer melt in a narrow 
channel. 

In the most general case, the system is composed 
of a blend of C distinguishable copolymer species in a 
fixed volume, V . Each copolymer of species a is itself 
constructed from distinct monomer types (such as 
polystyrene, polyethylene, etc.) Each copolymer species 
a is also taken to be monodisperse with a polymerization 
index Nq^ and an entanglement length Neo^. In addition, 
the system contains a set of rigid walls which are con- 
structed from distinguishable solid materials. 

We derive a hydrodynamic model for multiple vis- 
coelastic fluids at low Reynolds number by employing 
the principles of irreversible thermodynamics as outlined 
in 11]. The method described therein consists of a 
formal procedure for constructing a Rayleigh functional 
IZ = F + {1/2)W where F represents the time deriva- 
tive of the free energy and W is a dissipation functional. 
Extremizing IZ with respect to the component velocities 
produces hydrodynamic equations of motion while pre- 
serving the correct Onsager couplings. Due to space con- 
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straints, we present the results of the derivation here, 
deferring the details to a later publication. 

The fluid composition at a given point r is described 
by monomer volume fraction fields (paii^) and wall vol- 
ume fraction fields '^j(r). The first index on 
indicates the copolymer species, and the second index 
indicates the material type (polystyrene, polyethylene, 
etc.) Note that in general, polystyrene monomers as- 
sociated with one copolymer species ^ii(r) exhibit dif- 
ferent dynamics from polystyrene monomers associated 
with a second copolymer species ^21(1*) requiring us to 
treat them as distinct fiuids. The field '^j(r) indicates 
the volume fraction of solid material of type j at loca- 
tion r. By definition, all volume fraction fields sum to 
unity, T^ai^'ai-^^ji^j = 1- To save space, we employe the 

shorthand notation J^ai = — ^i=i- 

In the absence of chemical reactions, the number den- 
sity of each monomer species is conserved, as expressed 
by the continuity equations 



dt 



(1) 
(2) 

(pai and vj^(r) is the 



where v^^(r) is the velocity of fiuid 
velocity of the solid material 

Each velocity field may be split into two parts v^- = 

vt+w^.. The tube velocity vt = T^ai^ti^ti^T^j ^t^t 
represents the motion of the network of topological con- 
straints in an entangled polymer melt as described in 
Brochard's theory of mutual diffusion 12], and w^^ rep- 
resents the velocity of a given fiuid relative to vt- The 
stress division parameters Q^^i(r) and a'j (r) are obtained 
by balancing the frictional forces acting on the network 
as shown in glL] . 

The relative velocity fields are given by 



iVp- 



(3) 



which are produced as a result of imbalanced osmotic 
forces (/>Q;iV/iQ;i, pressure gradients Vp, viscoelastic forces 
V-cr, and external body forces fg. The relative velocity of 
each fiuid is inversely proprtional to its friction coefficient 
Cti = j^Coi^aiNa/Nea whcrc (oi IS the monomer friction 
coefficient of material i and vq is the volume occupied by 
a single monomer. 

The local force imposed on the ffuid by each solid ob- 
ject is 

f/(r) = + ^jVp - ap ■ a + (fivj - vt) (4) 

In the case where the wall velocities v'j are specified, 

we may integrate over f/(r) to obtain the net force and 
net torque on each object, potentially enabling numerical 
rheometric experiments. 



In the limit of low Reynolds number, the inert ial terms 
may be neglected, and the momentum transport equation 
becomes a force balance equation. 



= Vtt + Vp - V • o- - f ^ - f ^ 



(5) 



This equation implicitly determines the mean velocity 
field V — Tiaicl^ai^f + ^j^j^f which balances the net 

osmotic force Vtt = Tjo,i(t)ai^ 1-4 ^"^j^j^ 1^^ •> pressure 
gradient Vp, the viscoelastic forces V • tr, the net wall 
force f^, and the net body force acting on the fiuid . 

Energy is dissipated in the system by each of the terms 
in the dissipation functional W^^], including friction due 
to the motion of each fluid relative to the polymer net- 
work Ccxi{^ai — Vt), frictiou between the fluid and the 
walls Cfi'^f — Vt), and elastic deformation of the poly- 
mer network a : Vvt- Note that in the absence of vis- 
coelastic effects, the quantity V • cr reduces to the usual 
viscous dissipation term VjV'^v. 

The mean velocity v and pressure fields p must simul- 
taneously satisfy the force balance condition Eq. and 
the mass conservation condition, V • v = 0. These rela- 
tionships are not easily inverted and are obtained using 
an iterative numerical procedure. Once v and are 
known, the tube velocity may be calculated using the 
relationship 



Vt 
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Constitutive equations and are required to mea- 
sure the bulk elastic stresses cr^^ and shear elastic stresses 
cr^^ induced by deformations of each entangled polymeric 
component. 
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Although there are many possible constitutive equations, 
we have chosen to employ the Oldyroyd-B model in which 
the total viscoelastic force is produced by the sum of the 
elastic forces contributed by each polymeric component 
and a dissipative viscous force, V • tr = Xl^il^ ' + 
^ ' + r]sV^VT' This model allows us to study a 
purely viscous fluid {Gai = K^i = 0), a purely viscoelas- 
tic melt {r]s = 0), or anything in-between. The deriva- 
tives in {TJ and ^ are upper convected time derivatives, 

a - dtcr-\-VT ■ Vtr — a • Vvt — (Vvt)^ • cr. and the shear 
strain rate tensor is k!^ = diV^rj. + djv'^ ~ 1 ' 
The concentration dependent bulk moduli K^i and shear 
moduli Gai measure the stresses imposed on each compo- 
nent by the applied strains. Following Tanaka's example 
9], we employ moduli of the form Gai{<p) = Goai<Pai 
Koci{(l)) = KoaiO {(l)ai — fai) whcrc is the step func- 
tion. The quantity Tq, represents the reptative disen- 
tanglement time associated with the characteristic en- 
tanglement length A^e of copolymer a. For simulations 



with large elastic strains, this constitutive equation may 
be replaced with a more sophisticated phenomenological 
model or one based upon the SCFT microphysics 13]. 

A self consistent calculation produces the chemical po- 
tential gradients which drive phase separation. The free 
energy of the melt is F = U — TS where the enthalpy is 
given by 



U/kT: 
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(9) ■ 

The arid denote segment-segment and segment- 
wall interactions, respectively, and the entropy is 



(10) 



The quantity Qa is the partition function over all config- 
urations of a single copolymer species a in the presence 
of the fields co^-. 

Dynamic SCFT methods postulate that the system 
variables may be divided into two sets which relax on 
widely separated time scales. A "slow" set of parame- 
ters includes the elastic stresses and the conserved fields, 
mass, momentum, and monomer number. The remain- 
ing variables, including the conjugate chemical fields 
co'^- are "fast" parameters, which satisfy the local ther- 
modynamic equilibrium (LTE) conditions (f)ai[<^] = (j^ai 
T3]. which are derived by extremizing the free energy 
with respect to the conjugate fields uo'^^ir). The non- 
equilibrium volume fraction fields (pai and give rise 
to imbalanced chemical potentials /i^^(r) 
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which drive diffusion and 



phase separation. 

Calculations are performed in a periodic cell in two or 
three dimensions in which non-periodic geometries may 
be constructed using rigid materials. Spatial derivatives 
are computed pseudo-spectrally in Fourier space which 
enables the resolution of sharp domain interfaces with 
a small number of grid points and lends itself readily 
to parallelization. The multi-fiuid model is solved iter- 
atively for the cr,w, f^,v and p fields. Updated pres- 
sure and velocity fields are obtained from a projection 
method and frictional forces penalize fiow relative to the 
walls. The method repeats until a pseudo-steady state 
is achieved which satisfies the no-fiow, no-slip, continu- 
uity, and force balance conditions (|5|). The new veloc- 
ity fields are then used to transport the volume fraction 
fields. Similarly, the local thermodynamic equilibrium 
conditions are solved by the iterative application of a hill 
climbing technique. 

For concreteness, we have chosen to illustrate the tech- 
nique by examining a single representative system in de- 
tail, namely the phase separation of a quenched triblock 




FIG. 1: [top] Time evolution of a quenched triblock melt 
in a sib-micron channel at several flow rates as characterized 
by the dimensionless Weissenberg number, (a) Wi = 0.005 
(b) Wi = 0.013 (c) Wi = 0.029 (d) Wi = 0.108 [bottom] 
Lamellae orientation vs. time. The inset shows that time 
required to achieve Oxy — 0.40 diverges as the Weissenberg 
number approaches Wi = 0.02. 



copolymer melt in a sub-micron channel. The system 
consists of an incompressible melt of a symmetric ABC 
triblock copolymers (/i = 1/3) quenched from a disor- 
dered state into the intermediate segregation regime, in 
a channel of diameter D = 21 Rg. The Flory incompat- 
ibility parameters are fixed at x^ = 25, and = 12 
and all remaining parameters were chosen to approx- 
imate a polystyrene melt, with Co = 10~^Ns/m and 
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FIG. 2: Examples of 2D HSCFT systems: (a) A homopoly- 
mer blend A+B in two concentric, rotating cylinders (b) 
Nano-rods and nano-spheres in a phase separating blend (c) 
Viscoelastic flow of a blend through a microfluidic rectifier (d) 
An off critical diblock (/ = 0.3) in a micro-reactor channel. 

Go 2000kPa E^. The flow rate is characterized by 
the dimensionless Weissenberg number Wi = r|vmax|/^ 
where v^ax is the maximum mean velocity. 

The images in Fig. Q illustrate the morphology as a 
function of time at various Weissenberg numbers, where 
solid material is black and the copolymer blocks are red 
(A), green (B), and blue (C) respectively. We observe 
two distinct steady state configurations consisting of wall 
parallel and wall perpendicular lamellae, as illustrated 
by sequences (a) and (d) at t = 90. In the absence of 
flow, boundary wetting conditions dominate, resulting in 
a perpendicular orientation, and at high flow rates shear- 
ing effects dominate, resulting in a parallel orientation 
with blocks arranged in a CBA ABC configuration. The 
two effects compete, and we anticipate that there is some 
minimum Wiessenberg number Wic required to achieve 
a wall parallel configuration. 

To quantify these observations, we construct an orien- 
tation parameter o^^^y J dO s {0) [cos'^ O—sin^ 0]/{JdOs{0)) 
which measures the scattering difference between the x 
and y directions in Fourier space, where s{0) is the ra- 
dially averaged scattering intensity. The orientation pa- 
rameter is plotted as a function of time for several val- 
ues of Wi in Fig. ^ where Oxy = 1 corresponds to wall- 
parallel lamellae and Oxy = — 1 corresponds to the wall- 
perpendicular conformation. 

We estimate the value of the critical Weissenberg num- 
ber, Wic, by plotting the time needed to achieve a par- 
allel orientation, defined by Oxy{t) > 0.40, as illustrated 
by the inset in Fig. ^ It is clear from the figure that 
the alignment time diverges near the critical Weissenberg 
number, allowing us to estimate Wic ^ 0.02 in a chan- 
nel of diameter D = 21Rg. Furthermore, by comparing 
the two slowest cases, Wi = 0.005 and Wi = 0.013, we 
observe that a flowing melt with Wi <C Wic approaches 
its steady state configuration more rapidly than a quies- 
cent melt. Similarly, by comparing the Wi = 0.053 and 
Wi = 0.108 cases, we observe that a higher shear rate 
produces more rapid alignment in the Wi Wic regime. 

This numerical experiment also illustrates a type of 
dynamic asymmetry which arises in triblock copolymer 
melts. Although the system is compositionally symmet- 
ric, /a = Ib = fc^ we see distinctly asymmetric struc- 



tures in the early stages of phase separation as illustrated 
in FigH At t ~ 1, we find wall parallel structures pro- 
duced by preferential wetting of the walls by the A and 
C blocks, and at t ^ 5, we see an asymmetric drop- 
in-matrix morphology, indicating that early growth of 
density fluctuations is suppressed in regions which are 
enriched in block B segments. Since the enthalpic terms 
are all identical, we conclude that these effects are pro- 
duced solely due to differences in the conformational en- 
tropy between the copolymer blocks. To our knowledge 
this type of dynamic asymmetry has not previously been 
discussed. 

In addition to the example provided above, HSCFT is 
capable of modeling systems of far greater complexity. As 
illustrated in Fig. [2j some of the method's more unique 
capabilities include simulation of (a) the effects of moving 
machine parts on polymeric fluids (b) the unconstrained 
motion of solid nanoparticles (c) viscoelastic flow through 
irregular geometries and (d) flow in microfluidic channels. 

In summary, we have developed a mesoscale technique 
which models the effects of viscoelastic, Navier-Stokes 
type hydrodynamic transport on inhomogeneous poly- 
meric fluids in arbitrarily complex channels. In contrast 
with DSCFT methods, HSCFT is capable of describ- 
ing the influence of moving machine elements, nontrivial 
pressure induced and drag induced velocity fields, and the 
unconstrained motion of solid nano-particles. We believe 
this method has the potential for broad application and 
represents a first step toward more realistic simulations of 
complex fluid flows in industrial processing applications. 
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